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We consider a one-dimensional Josephson junction array, in the regime where the junction charg- 
ing energy is much greater than the charging energy of the superconducting islands. In this regime 
we critically reexamine the continuum limit description and establish the relationship between pa- 
rameters of the array and the ones of the resulting sine-Gordon model. The later model is formulated 
in terms of quasi-charge. We argue that despite arguments to the contrary in the literature, such 
quasi-charge sine-Gordon description remains valid in the vicinity of the phase transition between 
the insulating and the superconducting phases. We also discuss the effects of random background 
charges, which are always present in experimental realizations of such arrays. 
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I. INTRODUCTION 

The standard description of Josephson junction ar- 
rays isgiven by the theory by Bradley and Doniach 
(BD) [y. This theory predicts that a Josephson junc- 
tion array can either be in superconducting or insulating 
regime, depending on the ratio of the Josephson to the 
charging energy of the superconducting grains. The two 
phases are separated by the Kosterlitz-Thouless transi- 
tion. 

The theory of Bradley and Doniach assumes that the 
gate voltage applied to the arrays is equal to zero. This 
implies that the electrostatic energy of the islands is 
minimized when their charge is equal to the integer 
number of the elementary charge of Cooper pair. More 
recent work of Glazman and Larkin (GL) [2| demon- 
strates that the gate voltage could be another impor- 
tant parameter of the Josephson junction arrays, as two 
additional phases of the array are possible when it is 
varied. 

The BD transition was recently observed in one di- 
mensional Josephson junction arrays However, the 
results of Ref. [3| show only qualitative agreement with 
the theories of Refs. U and |3|. 

First of all, the BD theory completely neglects the 
capacitance of the Josephson junctions C as opposed 
to the capacitance of the superconducting islands Cq. 
GL does take nonzero C into account, but only as a 
small perturbation. In experiment, however, C 3> Cq. 
Not only does it mean that the junction capacitance 
terms are large, but it also leads to the long ranged 
Coulomb interactions between grains which are far away 
from each other, which is something both BD and GL 
neglect. As a result, the point in the parameter space 
where the transition happens is experimentally far from 
the theoretical point predicted by the BD theory. 

It is therefore of some interest to develop a theory 
of the Josephson junction arrays which would take the 
condition C 3> Co into account. In fact, such theories 
were already proposed in the context of both tunnel 
Q and Josephson [?| junction arrays. The main idea 
of this approach is to concentrate on the dynamics of 
charge of the islands as opposed to their superconduct- 
ing phases. As a starting point of the theory, one takes 



the problem of single Josephson junction with nonzero 
capacitance C treated in the full quantum mechanical 
fashion. The solution of this problem is a Josephson 
junction described by its quasicharge @. For a system 
of connected junctions this description in the continu- 
ous limit yields the sine- Gordon Lagrangian describing 
dynamics of the quasicharge. 

The sine-Gordon model in (l-fl)-dimcnsions is one of 
the best studied models of field theory. This model ad- 
mits a quantum phase transition: the excitation spec- 
trum depends on the coupling constant and becomes 
gapless when this constant exceeds a certain critical 
value. In the underlying Josephson array this transi- 
tion corresponds to the superconductor-insulator tran- 
sition predicted by BD. We have to note, however, that 
according to the opinion widely circulating in the lit- 
erature the quasicharge sine-Gordon description of the 
Josephson junction arrays loses its validity at large cou- 
pling constants and therefore cannot ascertain the ex- 
istence of the transition p| . 

In this paper we critically examine the conditions 
of validity of the sine-Gordon description of Joseph- 
son junction arrays. We conclude that the quasi-charge 
sine-Gordon description remains self-consistent in the 
region of large couplings. Hence one can use it to get 
quantitative information about the array even in the 
vicinity of the transition. 

We have to emphasize that the BD theory also de- 
scribes the transition in terms of the sine-Gordon equa- 
tion. However, being written in the limit C Co, the 
parameters of that equation are different from the one 
extracted from the sine-Gordon equation considered in 
this paper. We conclude therefore that while the univer- 
sality class of the transition does not depend on whether 
C is large or small, the details of the description do de- 
pend on it. 

Secondly, contrary to Ref. [2| , the experiments failed 
to observe any additional phases of the array due to the 
change in gate voltage. Most likely, this is the result 
of the presence of random background charges in the 
realistic Josephson junction arrays, which lead to the 
voltage being randomized along the length of the array 
0- An important outstanding question is thus the ef- 
fect the random charges may have on the array and how 
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their presence affects the transition. 

The problem of Josephson junction arrays with ran- 
dom background charges is equivalent to the problem 
of randomly pinned charge density waves Q , which was 
extensively studied in the literature. The best under- 
stood case is that of classical pinned charge density 
waves, which corresponds to the array deeply in the 
insulating phase. In that regime, the random back- 
ground charges lead to the AC conductance of the array 
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on the condition that the array length is much longer 
than the so-called Larkin length l p 12J. Here u> p is 
the pinning frequency. Additionally, for a special value 
of the sine-Gordon coupling parameter /3 2 = 4n (see 
later in the text for a precise definition of /?), the AC 
conductance is known to go as [Ll |: 



<j{oj) oc w 2 In w 



(2) 



This can in principle be checked experimentally. In this 
paper we discuss a relationship between 0, l p , oj p and 
parameters of the junction array. 

Unfortunately, at other values of the coupling param- 
eter (3, the solution to the random background charge 
problem is not known. Therefore, we are not able 
at this time to discuss how the presence of the ran- 
dom background charges affects the superconductor- 
insulator transition. 

The rest of this paper is organized as follows. 
In the next chapter we derive the quasicharge sine- 
Gordon equation. Then we discuss its applicability 
and show that it can indeed describe the transition 
superconductor-insulator transition. In the end we dis- 
cuss how the random background charges affect the 
properties of the array. 



II. DERIVATION OF THE MODEL 

Let us consider a one-dimensional array of weakly 
coupled Josephson junctions (see Fig. 1). A single junc- 
tion is described by the Hamiltonian 
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is the charging energy of the junction, 



Ej is the Josephson energy and cj) is the phase difference 
on the junction. 

We work in the limit Ej 3> E c . In this regime, the 
energy spectrum consists of narrow bands separated by 
gaps (see Ref. 0,0). The interband splitting frequency 
is estimated as 



W = 2^E~ 



(4) 



We shall restrict our consideration to the lowest band. 
Therefore W will serve as the ultraviolet cut-off in our 
effective theory. The energy eigenvalues in the lowest 
band are labelled by the quasicharge q: 



E = E b [l- cos(27rg)] 
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FIG. 1: An equivalent scheme of an elementary cell of the 
array. Black dots denote the superconducting islands. 



Here 
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The junction array consists of junctions separated by 
superconducting islands. Following we treat the ar- 
ray in the adiabatic approximation, so that the n-th 
junction is characterized by a slowly varying function 
q n {t). The characteristic frequencies are supposed to 
be much smaller than the interband splitting W . The 
charge on the n-th junction is 2e x q n . Since the to- 
tal charge in the unit cell is zero, the charge on the 
wire connecting the two junctions is 2e(q n — q n +i)- The 
Coulomb energy of the wire is 



^Coulomb - -g^T- ~ qn+1 ^ ~ / 
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where a is the size of the island. The inductive energy 
of the junction is 



^Inductive - lj2^ 2e ^ L ^ 



(8) 



where L is the inductance of the islands. The exper- 
iments conducted in Ref. 0] indicate that real ar- 
rays have a very considerable inductance (L 3> J? E ) , 
though its microscopic origin remains somewhat ob- 
scure. 

Combined, Ecis. l5l7l8fl produce the Lagrangian de- 
scribing the array, 



£ = dx 
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By introducing a new variable Q — 2nq/P the latter 
Lagrangian can be rewritten in the form of the sine- 
Gordon model: 



fiv r . 



dx 



1 



where 



ml 



7T 2 E b 

e 2 L ' 



m| (1 

f3 2 y 



cos j3Q) 

(10) 

(11) 
(12) 



3 



and 



/3 2 = 
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where Rq = h/4e 2 is the quantum of Cooper pair re- 
sistance. 

The sine-Gordon equation (|10fl goes through a quan- 
tum phase transition as /3 2 is tuned through 87r. Let 
us estimate whether this value of /3 can realistically be 
achieved. Assuming that the island has a form of a 
th in sup erconducting wire we can use for its impedance 
\J LjC^ the formula derived by Likharev : 
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where A is the penetration depth, A is the cross-section 
area of the wire, d is its thickness and D is the distance 
from the gate. Taking into account that hc/e 2 « 137 
we get 
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We see that /3 2 can be driven through the transition by 
changing A which can be achieved by adding nonmag- 
netic impurities |17|. 



III. APPLICABILITY LIMITS OF THE 
SINE-GORDON MODEL 

As we have mentioned above, the sine- Gordon model 
Eq. (|lf)|> is one of the best studied models of one- 
dimensional field theory. Its spectrum and behavior of 
correlation functions dramatically depend on the value 
of coupling constant /3 2 , because 1 controls the renor- 
malization group (RG) dimension of the cosine term 
in Eq. I|1U|) . At 1 < 8ir, the cosine term is relevant, 
all excitations have spectral gaps and the system is an 
insulator. At f3 2 < 4tt the spectrum consists of electri- 
cally charged solitons and antisolitons and their neutral 
bound states, while at 8n > (3 2 > An the bound states 
disappear. At 1 > 8ir, the cosine term is irrelevant and 
the spectrum becomes gapless. This is a superconduct- 
ing regime; in the absence of disorder it allows ballistic 
transport through the system. 

Let us critically examine conditions of validity of the- 
ory UlUfl . It is based on the following assumptions: 

• The harmonicity of the effective potential. The 
potential acquires cosine form JSJ in the limit 
Ej 3> E c . The presence of higher harmonics will 
affect the integrability of the sine-Gordon model, 
but will not change the fact of the transition 
since these harmonics will become irrelevant even 
sooner than the primary one. 

• The abiabaticity. All characteristic energies of 
sine-Gordon model (in particular, the spectral 
gaps) must be much smaller than the interband 
splitting W of Eq. |QJ. 



• The continuous approximation. The discrete ar- 
ray is replaced by the continuous one. This ap- 
proximation requires that characteristic wave vec- 
tors are much smaller than 1/a and the energy 
gaps are much smaller than A = hv c /a. 

• Absence of dissipation. It is assumed that islands 
are superconducting and there are no normal re- 
sistors in the scheme. 

• We assume that all characteristic frequencies are 
much smaller than the plasma frequency in the 
islands so that electric charge spreads instanta- 
neously through the entire island. 

Let us look closer at the requirement of adiabatic- 
ity and validity of the continuous approximation. The 
spectrum of the sine-Gordon model consists of particle 
branches with the relativistic dispertion 
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where the spectrum consists of solitons s and antisoli- 
tons s with m s — iris — M and (at /3 2 < 47r) their 
bound states (breathers) with spectral gaps 
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The largest spectral gap is of order of M. The self- 
consistency of the theory requires that its value must 
not exceed the cut-off. According to 16j, we have 
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16A 2 r(l + /3 2 /8tt) 



M _ 

T ~ 0Fr(l/2 + 7 /2) 
The following inequalities must be fulfilled: 
M < A < W 
Inequality M <C A is equivalent either to 
too < 7 A, E b < l/e 2 L 
at small /3 2 or to 

too < A, E b < e 2 /Cb 



2-,3 2 /47i 



(18) 



(19) 



(20) 



(21) 



at (3 2 « 8tt. The validity of the continuous approxima- 
tion requires that A <§; W, that is 



< EjE e , or (3 2 < WCo/Ae 2 



(22) 



From the previous inequalities it follows that the 

only condition on the latter quantity is WCo/e 2 <C 
i 

■g 2 -^ 4 exp[2^/Ej / E c ] which still leaves room for /3 2 pa 



IV. PHASE TRANSITION 

As we have mentioned above, one way to drive the 
system through the transition is by changing f3 2 . In 
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fact, it is also possible to change or equivalently mo, 
to go through the transition. This would correspond 
to the experimental set up Q, where Ej is lowered by 
applying the external magnetic field to the junctions in 
the SQUID geometry, leading to the increase of E^ and 
mo. 

The Renormalization Group (RG) flow diagram of 
the sine-Gordon equation is shown on Fig. 2. If 
/3 2 < 8tt, the sine-Gordon equation is always in the 
massive regime. This corresponds to the insulating be- 
havior of the array. However if (3 2 > 8ir, two regimes are 
possible. For sufficiently large mo, the behavior is still 
massive. For mo smaller than a certain critical value 
m c cx A (/3 2 — 8n) , the cosine term in <(TT)|) becomes ir- 
relevant and can be neglected. The sine-Gordon equa- 
tion becomes massless. This corresponds to the part of 
the diagram whose RG flow lines end on the mo = 
axis. In this regime the charge density wave propagates 
ballistically. 



massive regime 




massless regime 



3 8ji 



FIG. 2: The RG flow diagram of the sine-Gordon equation. 
One can tune through the transition by either changing the 
value of (3 or mo. 

One of the more experimentally relevant predictions 
of the theory for the ballistic regime follows from ne- 
glecting the cosine term in Eq. I|I(J|I . Then the Joseph- 
son junction array can be described by the following 
very simple Lagrangian 



C = W c 



dx 
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This is none other than the so-called Luttinger liquid 
with the Luttinger parameter g = j^. It is well known 
[l8| that its conductance is given by 
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As (3 is lowered (or mo is increased, according to Fig 2), 
the conductance is decreased until it reaches its critical 
value (at (3 2 = 8tt), 
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If (3 is decreased further, the array enters the regime of 
the Coulomb blockade and becomes an insulator. 



V. RANDOM BACKGROUND CHARGES 

In the presence of random background charges, the 
charging energy of the superconducting islands Eq. 
gets modified, 



^Coulomb = J dx a 



^(q x -V(x)f, (26) 
zOo 



where V(x) is a time independent random function of 
the coordinate x with short ranged correlations 



{V{x)V(y)) = V 6(x-x'). 



(27) 



It is convenient to shift the variables q(x) — > q(x) + 
f Q dy V(y) to find the following Lagrangian for the 
Josephson junction array, 



C =2 
where 



dx 



-Q 2 
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X (x) = 2tt / dy V(y). 
Jo 



(28) 
(29) 



It is clear that x( x ) is a Brownian motion; \/ (x 2 ( x ) ) = 
2ny/xVo. At distances of the order of I cx V the 



phases of the cosine term in Eq. (|28|l become com- 
pletely uncorrelated, so at length scale much bigger 
than that we can consider cos{x(a;)} to be white noise 
in space. / is the length at which the random back- 
ground charge accumulates to be of the order of the 
Cooper pair charge. 

Under these conditions, the problem defined in 
Eq. H28fl becomes equivalent to the pinned charge den- 
sity wave problem studied in a number of publications 
[Uj • A particularly well understood regime is that of a 
Coulomb blockade, (3 <C 87r. There the sine-Gordon La- 
grangian Eq. I|28|l becomes purely classical, correspond- 
ing to the sine-Gordon equation 

2 

v -id?Q-v c d 2 x Q-^ S m{Q + X } = 0. (30) 

P 

If Qo(x) is a time independent solution to Eq. ij^Hjl . then 
the small oscillations around that solution are described 

by 

v c d 2 JQ + ml cos {Q + x} 5Q = v~ x lo 2 8Q, (31) 

where SQ is the amplitude of the oscillations and lu is 
their frequency. 

It is well known in the literature 0, IToL [TT| that this 
problem possesses a fundamental frequency called the 
pinning frequency u> p . At lu 3> oj p , it is possible to 
neglect the cosine term in Eq. (|31|l to find that the os- 
cillations SQ are the plain waves with the wave vector 
k and with the frequencies u> oc v c k. At lu <C lu p , the os- 
cillations are localized in space, and their spacial extent 
(localization length) l p is called the Larkin (or pinning) 
length. The pinning frequency is known to scale with 
mo as 
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(32) 
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while the Larkin length scales as 

1 -i 
lp = — ex m 3 . (33) 

LUp 

On the condition that the Larkin length is much smaller 
than the total length of the array, the conductivity of 
the array at a frequency to much smaller than the pin- 
ning frequency is given by 

<j(w) (X p{(jj), p{iS) OC Ld 4 , LU <C UJp. (34) 

Here p{u)) oc cj 4 is the probability that there exists a 
solution to Eq. (|31|l at a frequency uj (or the density 
of states of Eq. if3*TJ) h The formula Eq. J2H as well as 
Eq. I1M2I) could be checked experimentally. 



Finally, at a special value of 1 = Air, the Eq. (|28H can 
be solved using a completely different technique [l3j , to 
lead to the AC conductivity Eq. J2J. 

At values of f3 other than /3 2 <C &tt or (3 2 = An, the 
problem defined in Eq. I|28|l remains unsolved. 
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